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The current-induced resonant excitation of a magnetic vortex core is investigated by means of an- 
alytical and micromagnetic calculations. We find that the radius and the phase shift of the resonant 
motion are not correctly described by the analytical equations because of the dynamic distortion of 
a vortex core. In contrast, the initial tilting angle of a vortex core is free from the distortion and 
determined by the nonadiabaticity of the spin torque. It is insensitive to experimentally uncontrol- 
lable current-induced in-plane Oersted field. We propose that a time-resolved imaging of the very 
initial trajectory of a core is essential to experimentally estimate the nonadiabaticity. 



A spin-polarized current can exert torque to a fer- 
romagnet by transferring spin-angular momentum, i.e. 
spin-transfer torque. The spin-transfer torque provides 
full magnetization reversal, steady-state precession mo- 
tion, and domain wall movement [H, Q ■ It is composed of 
adiabatic and nonadiabatic spin torque terms in continu- 
ously varying magnetization. The adiabatic spin torque 
arises from the conduction electron spin whose projection 
on the film plane follows the direction of a local magne- 
tization, whereas the nonadiabatic torque arises from a 
mismatch of the direction as a result of the momentum 
transfer or the spin relaxation [1, 0, ■ 

Until now, the experimental threshold current density 
Jc to steadily move a domain wall has been reported 
to be about 10* A/cm^, too large for an application. In 
addition to the resonant depinning Q and the use of 
perpendicular anisotropy Q, an increase of (3, the ra- 
tio of the nonadiabatic spin torque to the adiabatic one 
can reduce Jc [I]- Despite its importance, however, 
the exact value of 5 is still under debate even in theo- 
ries [i, 0, [nl, H [ll; /3 = 0, /3 = a, and /3 ^ a where 
a is the Gilbert damping constant. This debate is also 
related to which mechanism between the Landau-Lifshitz 
damping and the Gilbert one is correct to describe the en- 
ergy dissipation 1J| . Experimental determination of /3 is 
essential to resolve this debate, but experimentally esti- 
mated values are also distributed; /3 = 8a Q , /3 = a [lB| , 
(}> a [III, (} = 2a [i3|,and /3 7^ a [11. Since most ex- 
periments have used the same material (Permalloy) , this 
wide distribution is caused by origins irrelevant to the 
material itself. 

The wide distribution can originate from the Joule 
heating, the edge roughness of nanowire, and the in- 
plane component of the current-induced Oersted field 
i/g". The Joule heating significantly affects Jc and wall 
velocity [H, IH. Therefore, it is difficult to precisely es- 
timate (3 when the Joule heating is not negligible, i.e. 
Jc > IQ^A/crn? . In a magnetic nanowire, the edge 
roughness distorts the domain wall [lO] and thus pre- 
vents a proper interpretation of experimental data us- 



ing theories derived for an ideal nanowire. A way to 
avoid the above issues is to experimentally study resonant 
motions of a magnetic vortex core (VC) in a patterned 
disk by injecting an alternating current of the order of 
lO^A/cm^. The magnetic vortex is an ideal system for 
the resonant motion study since VC can be considered as 
a topological point charge which efficiently responses to 
external forces 
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|22| . |23| | . It was experimentally con- 
firmed that the VC can be resonantly excited by an ac 
current [iHHIii. Even in this case, however, a very 
small in-plane component of the ac current-induced alter- 
nating Oersted field i/^" inhibits a precise estimation of 



the spin torque parameters [25|, l26| . Note that the Hq is 



not a current-induced field along the thickness direction 
of the disk, but an in-plane field caused by any geometri- 
cal symmetry-breaking of the system. The driving force 
due to the iJ^" of only 0.3Oe is as large as 30% of the to- 
tal resonant excitation ^2F| . Such a small iJ^" is difficult 
to remove since it is caused by an uncontrollable nonuni- 
form current distribution due to a geometrical symmetry- 
breaking such as electric contacts or notches. Therefore, 
an experimental way to estimate the (3 which is safe from 
the Joule heating and the edge roughness, and also in- 
sensitive to the HqI is highly desired. 

In this letter, we propose that a direct imaging of the 
very initial trajectory of VC induced by an ac current is a 
plausible way to experimentally estimate j3. We find that 
(3 does not change the resonant frequency, but affects the 
phase of resonant motion. The phase-shift is /3-dependent 
since (3 determines the tilting angle of very initial trajec- 
tory measured from the direction of the electron-flow. 
On the other hand, the HqI with a typical magnitude 
does not change the tilting angle although it affects the 
steady resonant motion. More importantly, the initial 
tilting angle is only one physical observable which can be 
directly compared to the analytical result, whereas the 
others such as the radius and the phase shift are not cor- 
rectly described by the analytical equations because of 
the dynamic distortion of VC. 

The current-induced motion of VC is calculated using 
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FIG. 1: (Color online) (a) Schematics of the model system, 
(b) Magnetization of a vortex in its initial equilibrium state. 
The height denotes the z component, whereas the color scale 
corresponds to the direction of the x-component of magneti- 
zation. 



the Thiele's equation with the spin-transfer torque terms 

(Eq. (HI)) [i,^,!!!,!!^. 



G(p) X (u - r) 



SU{r) 
Sr 



aDr + (3Du (1) 



where G{p) — —GopSz is the polarity (p± 1) dependent 
gyrovector, Gq is obtained from the spin texture as 



Go — 



Ms 



dVsm{0){V9 X Vip)-e^, 



(2) 



6 (ip) is the out-of-plane (in-plane) angle of the magneti- 
zation, Ms is the saturation magnetization, 7 is the gy- 
romagnetic ratio, u = uoexp{iujt)ey, uo{— ^bJP/gMs) 
is the amplitude of adiabatic spin torque, uj is the angu- 
lar frequency of the ac current, r{t) = X{t) + Y{t) is the 
time-dependent position vector of VC, and U{r) is the 
potential well. The damping tensor D is also obtained 
from the integration of spin texture as 



D = -■ 



M. 



'-is 
7 



dV[{V9V9 + sm'^{9)Vi^Vi^)]. (3) 



When VC is at the static equilibrium position, Gq is 
2iTLMsh and D is Go ln(i?/5)/2 where L is the thick- 
ness of disk, R is the vortex radius, and 5 is the core 
diameter. 

With X{t) = Xq exp(iw<), Y{t) = Yq exp{iLut) and 1 + 
{aD/G[)Y ~ Ij the solutions are in the form of X{t) = 
Xx cos(wt)— X2 sin(a;t) and Y{t) = Yi cos(a;t)— y2 sin(a;i) 
where 



Xi 
X2 
Yi 
Y2 



AujG[{I3 - a){ujl ~ uj^) + 2auj% 
AuJrC[{Lul - uo^) + 2a(l + C^a[3)Lu^], 
Au[{l + C'^a(i){ujl - t^2) - 2C'^a[3uj'^^]. (4) 



Here, A = uq/[{uoI — u?)'^ + (2GaWrw)^], = k/Gq is 
the resonance frequency, k = (dU / dr) jr is the effective 



stiffness coefficient of the potential well, and C \sD/Gq = 
\n{R/5)/2. From the eqs. dH), one finds the radius a{t) = 
\J X{tY + Y{tY, and the phase shift (j) between the phase 
of the core gyration and that of the ac current; 



— tan 



= tan 



1 - {iOr/uj f + 2C^a{p ~ a) 



G/3((wJw)2 -1) + 2Ga(l + C^al3) 



.(5) 



To verify the validity of the analytical solutions, the 
micromagnetic simulation is also performed by means of 
the Landau-Lifshitz-Gilbert equation including the spin 
torque terms; 



dM „ a dM 

— ^ - ^^xHe// + ^Mx — 

-I- woVM - f^M X VM 
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where ^eff is the effective field including the exter- 
nal, the magnetostatic, the exchange, and the current- 
induced Oersted field. The model system is a circular 
Permalloy disk with the thickness of lOnm and the di- 
ameter of 270nm (Fig. 1 (a)). The unit cell size is 
2 X 2 X lOnm^. The ac current with the maximum cur- 
rent density of 1.25 x lO"^ A/cm^ is applied along the y- 
axis. Standard material parameters for Permalloy are 
used; Ms = SOOemu/cm^, 7 = 1.76 x 10^Oe"is~\ 
a — 0.01, P — 0.7, and the exchange constant A^x = 
1.3 X lO^^erg/cm. 

First, we assume HqI = in order to investigate the 
effect originating exclusively from P on the resonant mo- 
tion. We will recall the effect of HqI in the last part. 
In Fig. 2 (a) and (b), we show analytical and modeling 
results of the time averaged value {Y{t) ■ Jq sin(wi))//o 
at various /3-terms. To obtain analytical results, we use 
G = 1.3 because R is 135nTO and S is lOnm, determined 
from the micromagnetic configuration at the initial equi- 
librium state. {Y{t) ■ Iq sm{ujt)) / Iq shows a peak at the 
resonance frequency Ur of 360MHz. ojr does not change 
with (3 whereas the peak structure becomes more asym- 
metric as P increases. In spite of qualitative agreement, 
however, analytical results are quantitatively different 
from modeling ones. This is because the radius of gy- 
roscopic motion is different between the modeling result 
and the analytic solution (inset of Fig. 2(b)). We at- 
tribute this difference to a dynamic distortion of VC. As 
shown in Fig. 2(c), the VC shape in the initial equilib- 
rium state is symmetric whereas it in a dynamic motion 
is asymmetric. The distortion changes the gyroscopic pa- 
rameter Go, the damping tensor D (thus, the parameter 
G) and the effective stiffness coefficient k since all pa- 
rameters are determined by details of the spin texture 
(See, eqs. ([2]), and the definition of k). From mi- 
cromagnetic spin textures of the vortex in the dynamic 
motion, we find that both G and k increase from the 
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FIG. 3: (Color online) (a) Phase shift as a function of the 
frequency for various /3-terms. (b) Initial vortex core trajec- 
tories for various /3-terms {uj — LJr)- Inset of (a) shows the 
vertical offset of the phase shift for various /3-terms. In the 
inset, open symbols are obtained from micromagnetic model- 
ing. Inset of (b) describes the definition of the initial tilting 
angle. 



FIG. 2: (Color online) The {Y{t) ■ lo sm{Lot)) / lo as a func- 
tion of the frequency obtained from (a) analytic solution 
of Thiele's equation and (b) micromagnetic simulation, (c) 
Comparison of the shape of vortex core, and (d) variation of 
Hi clS cL function of the gyration radius. The inset of (b) shows 
the vortex core trajectories. The inset of (d) shows variation 
of the parameter C (= D/GO) with the time evolution. Both 
insets are obtained ai P — and u — Ur- 



initial equilibrium values because of the dynamic distor- 
tion (Fig. 2(d)). The increase of k is much larger than 
that of C, and responsible for the reduced radius in the 
modeling results. This increase of k occurs at the very 
initial time stage, indicating that the assumption of the 
rigid VC and potential well is invalid except for the very 
initial trajectory. Therefore, even when -/Jq" is zero, it 
may be difficult to deduce important parameters such as 
P by directly comparing analytical solutions with exper- 
imental measurements of the steady gyroscopic motion. 

Nevertheless, it is worthwhile investigating how /3 in- 
duces the asymmetry in the peak structure. Fig. 3(a) 
shows analytical phase shifts at various /?-terms. As p in- 
creases, the absolute value of the phase shift (j) decreases 
(increases) for the frequency smaller (higher) than ujr- 
In other words, a vertical offset 0(/3) — 0(/3 = a) of the 
phase shift increases with increasing /3. From the eq. ^ 
with a <C 1 and /3 <C 1, one can find that the vertical 
offset is approximately tan~^[C(/3 — a)] and thus depen- 
dent on f3. This is why the peak structure becomes more 
asymmetric as increases. However, quantitative disagree- 
ment between analytic solution and modeling result was 
again observed (inset of Fig. 3(a)). The difference be- 
comes larger as the radius of core gyration increases, i.e. 
the frequency approaches LOr- This is also caused by the 
dynamic distortion as explained above. 



The vertical offset is /3-dependent since the initial tilt- 
ing angle Bint is determined by /3 (Fig. 3(b)). VC initially 
moves along the direction of the electron-flow. Because 
of the imbalance of magnetostatic field, VC experiences 
the centripetal force and starts to undergo a gyration mo- 
tion. In the absence of HqI, 9int of the initial trajectory 
can be obtained from the eq. ([T]) by dropping the po- 
tential gradient term since VC is initially at the bottom 
of the potential well where the gradient is zero. When 
HqI is nonzero, the potential gradient is no longer zero 
and could affect 9int- In order to investigate the effect 
of on Oint, we perform micromagnetic simulations 
for the initial trajectory with and without taking into 
account HqI. We assume that an alternating HqI is ap- 
plied along the x-axis and its magnitude is 0.3Oe which 
is similar with the estimated value in the experiment of 
the Ref. As shown in Fig. 4(a), the effect of H^^^ 

on the very initial trajectory is negligible whereas the 
difference in trajectories between the two cases becomes 
larger and larger as the time evolves. This insensitivity 
of the initial trajectory and thus Oint to HqI is valid for 
a different (3 (not shown). Thus, we drop the potential 
gradient term in the eq. ([1]) to derive Oint- With p = +1 
and the direction of the initial current along the -|-y-axis, 
one finds 



aDX-GoY = Gouo, 
GoX + aDY = -PDuo 



(7) 



where X and Y are the velocity along the x- and y-axis, 
respectively. By solving the eqs. ([7]) for X and Y, 



X = GqUqD- 



^2132 + Q2 ' 



Y = -uo 



Gl + a/3£)2 
a2D2 + d 



(8) 
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FIG. 4: (Color online) (a) Effect of in-plane Oersted field 
(^Oe) on the initial trajectory of vortex core (/3 = and 
uj = LUr), and (b) effect of /3 and C on the initial tilting angle. 



The initial tilting angle 0j„t is given by 



= tan"i (^) ^ (l 



C{ a ~ 13) 



(9) 



Note that the eq. ^ is equivalent to the equation for 
the vertical phase-shift with considering the sign of the 
initial ac current and a/3 <C 1 . It confirms that the verti- 
cal shift originates from the /3-dependent 9int- Fig. 4(b) 
shows 9int as a function of 13/ a for various values of C. It 
should be noted that Omtn obtained from modeling are 
in excellent agreement with analytical ones in contrast 
to the radius and the phase shift. It is because VC re- 
tains its equilibrium shape at the very initial time stage 
where C and k hardly changes. For the tested Permalloy 
disk (C = 1.3), the difference in Oint between = and 
(3 — IS about 6 degree which may be small to exper- 
imentally measure. However, the /3-dependent Oint be- 
comes larger as the parameter C increases. For instance, 
the difference in Oint between (3 — Q and /3 = So; is about 
18 degree at C = 4. The C increases as the disk diameter 
(core diameter) increases (decreases) . The core diameter 
5 decreases with decreasing the disk thickness [s^. The 
equilibrium value of C can be determined by micromag- 
netic calculation or MFM imaging for a given disk ge- 
ometry. Consequently, a time-resolved magnetic ima ging 
with high spatial resolution such as X-ray microscopy |3l| 
for observation of the very initial trajectory of VC in a 
wider and thinner disk is a possible way to experimen- 
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